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Recently dispersion relations have been applied to hard exclusive processes such as deeply virtual 
Compton scattering, and a holographic principle was proposed that maps out the generalized parton 
distributions entering the soft matrix elements for the processes from their values on a given kine- 
matical ridge. We examine possible pitfalls associated with the implicit, direct identification in this 
approach of the physical hadronic states with colored partons, and suggest an improved treatment 
of this assumption. 



A number of studies have recently advocated using 
Dispersion Relations (DRs) both to facilitate the extrac- 
tion from deeply virtual exclusive experiments, such as 
Deeply Virtual Compton Scattering (DVCS), of the soft 
matrix elements for hard exclusive processes, the Gen- 
eralized Parton Distributions (GPDs), and to determine 
their model parameters [TJ [51 [3] . In this paper we are go- 
ing to show that there are important limitations to the 
use of DRs for processes described by GPDs. 

DRs have a long history in hadronic physics. For a 
general exclusive, two body hadronic reaction, invariant 
amplitudes can be written in terms of energy and angle 
variables, such as the Mandelstam variables s and t, or 
V = {s — u) /AM and t. When the energy variable is con- 
tinued into the complex plane, the amplitudes become 
holomorphic functions, i.e. analytic functions over re- 
gions of the complex plane. Unitarity of the amplitudes 
determines the physical intermediate states that, in turn 
determine branch cuts in the complex energy plane. Each 
physical state has a kinematic threshold that fixes the 
branch point. 

DRs were derived for inclusive Deep Inelastic Scatter- 
ing (DIS) as well, when viewed as forward virtual Comp- 
ton scattering [11 . in this case translates into the virtual 
photon's energy in the laboratory system {v = VLab), 
in turn connected to the variable u = l/xsj, where 
XBj — /2Mi>. DIS can be considered as a special case 
of elastic scattering where unitarity relates the imagi- 
nary part of the forward amplitude to the total cross 
section, the inclusive sum over all physical final states 
allowed by the energy. As discussed thus far, DRs do 
not necessarily include the partonic structure of the tar- 
get. Partonic degrees of freedom are integrated over. In 
fact, all remaining kinematical variables, including XBj, 
can be considered to be fixed by the kinematic conditions 



'Electronic address: gary.goldstein@ tjifts.edu| 
tElectronic address: sl4y@virginia.edu 



"external" to any partonic loop or QCD elaboration. 

A connection with the partonic structure, through the 
Operator Product Expansion, therefore QCD, can be es- 
tablished e.g. by following the derivation in Ref. 0], 
where the important assumption is made that the sup- 
port for both the integrals defining the Mellin moments 
of the operators and the final amplitude is XBj G [—1, 1], 
in the asymptotic limit, ^ oo. It is also assumed that 
the intermediate states that are summed over in a fac- 
torized amplitude are physical. This leads to the identi- 
fication of the (twist two) quark distribution, HixBj) |1] 
with the measured structure function, F2{xBj), or the 
imaginary part of the forward amplitude. More specifi- 
cally, two steps are taken in establishing DRs for DIS: 
i) the identification of the physical threshold for the 
scattering process, z^f/j, with the continuum threshold, 
vc ~ M + rriTr, or xc — 1, appearing in the integral defi- 
nition of the scattering amplitude; ii) the identification of 
XBj with the partonic variable present in the factorized 
amplitude. 

In this Letter we argue that these assumptions cannot 
be carried straightforwardly to the off-forward case de- 
scribed e.g. in DVCS. In fact, as explained later, one has 
a mismatch between the supports for the scattering am- 
plitude and for the corresponding DR, namely vc 7^ vth- 
This mismatch is a straightforward consequence of t- 
dependent physical thresholds, not present in the DIS 
forward/elastic case, that are long known to hinder the 
useful and practical applications of DRs. The mismatch 
exists for both v and C,. This point does not touch upon 
the partonic aspects of the process. However, in the fac- 
torized form of DVCS, described by a handbag picture, 
XBj in ii) is replaced by two longitudinal fractions, X 
and C, where X = (kq)/{Pq) and the external variable, 
the skewness, C = {q^)/{qP) ~ Q"^ / {2MvLab), A being 
the momentum transfer for the two body scattering pro- 
cess, — t (we will use either the set {X,(^,t) or the 
alternative variables, {x = '^_^/2 ; ^ — 2^ ' ^) through- 
out the paper- see Ref. [HI [5] for reviews on DVCS). The 
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expression for the DVCS ampHtude in QCD factorization 
is 

flavors 

where the analog of the Compton Form Factor (CFF) is 



Hf{x,^,t) 
X — + ie 



(2) 



The GPD H{x,£^,t) is convoluted with the hard part, 
1 /{x — ^ + ie), and integrated over x in the range [—1, 1]. 
Crossing symmetry is implemented by 

Hf\x, ^, t) = Hf{x, t) T Hf{-x, t), (3) 

recalling that for PDFs, q{—x) = —q{x) relates negative 
X to positive x antiquark probability. 

It follows straightforwardly from Eq.([2]) that 
Im7i(^, t) = t). To relate this to the discontinuity 

across the physical branch cut of a holomorphic function, 
unitarity is invoked through the insertion of a complete 
set of intermediate states. 



ImH(C, t) ^ J dX [Six - + S{X)] (4) 

N 

X {2tt)S{XP+ +P+-P+) 



The resulting analytic structure allows the DR to be writ- 
ten, 



P.V. / dx 



a; - C 



(5) 



where the subscript unphys emphasizes that the integra- 
tion should be over the whole range, but because the 
integration variable is now interpreted as the skewness, 
"external" to the quark loop, a threshold mismatch en- 
sues due to the inelasticity of the two body process for 
non-zero t. In fact Qh [-t + {t^ - AAPtY/'^]/2M 
for » t, the physical threshold for the two body, 
1*P ~^ iv' scattering process, originates from the limiting 
values for the angles defining the invariant t = (q — q')^. 
One obtains in the limit » t, Unin = /^s—{q^^^ ~ 
q'CMy ^ _M'^(^/{l-i^). Notice that for DIS, the phys- 
ical and continuum thresholds coincide because the final 



photon has the same as the initial one, tmin = and 

Cth = Xth = 1 = Cc- 

In DVCS the region x e [^f^,, 1] is unphysical and the 
second term in Eq.(|5| cannot be obtained from experi- 
ment. The physical meaning of this discrepancy is illus- 
trated in Figjl] where both the continuum and physical 
thresholds for several variable describing DVCS, s, v and 
C, are plotted as a function of t. For s, as increases, 
only higher and higher invariant mass states are sam- 
pleded. Under s u crossing, there are corresponding 
branch cuts for negative v. So it is not clear how the dis- 
persion integral can be written in the partonic variables. 
Although the mismatch between physical and continuum 
thresholds addresses the issue of the physical interpre- 
tation of GPDs, it was a well known problem for two 
body scattering processes [7| where it was dealt with by 
either constructing models for the analytic continuation, 
or developing some other prescription. The threshold 
mismatch seen in these fixed t DRs in C, could be reduced 
by introducing new variables, a method used in hadronic 
processes. We will show the consequences of introducing 
a jet mass in the factorized picture [S]. 

To illustrate the different physics involved in forward 
and off-forward processes respectively, we discuss the 
proof of the DR given in Ref. [1 . This was obtained sim- 
ilarly to the DIS case (see Jaffe |1] for example). The 
hadronic tensor related to the forward Compton ampli- 
tude T^^'^{x,Q^) can be given a partonic interpretation 
when the operators are expressed via interacting quark 
fields and subjected to the OPE. The coefficients of the 
leading twist terms in that expansion are the Mellin mo- 
ments of the quark distribution functions H{x). Sum- 
ming this geometric series (for \x\ > 1) leads to the form 
of the DR for T'^'ix, Q^) with integrand H{a)/{x ± a). 
Because the Compton form factor is known to satisfy an- 
alyticity as a forward elastic amplitude, that analyticity 
allows the continuation of the integration to the com- 
plex X plane and the DR follows. In Ref. [T] the GPD, 
H(x,S,,t), enters observables through integration over x 
as in Eq.([2]). The denominator in the integrand, which 
arises from the light cone limit of the struck quark's prop- 
agator, can be written as a geometric series in |. Because 
the GPD must satisfy polynomiality in (the x moments 
are polynomials in ^ with i-dependent coefficients) , based 
on the underlying covariance, the resulting series must 
converge for large |^| > 1. So in the complex ^ plane the 
Ti.f{£,,t) will be analytic for the unphysical |^| > 1. 

Then, by analogy with the hadronic DR, it is assumed 
that there is a physical branch cut from -1 to -1-1 on the 
real ^ axis and no other poles or cuts. For this interpreta- 
tion however, the intermediate states, the s-channel cuts, 
have to be determined, given non-zero t and Q2 But 
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FIG. 1: (color online) Thresholds for the variables v and s 
plotted vs. t. The dashed lines are the continuum thresholds 
whereas the full lines represent the physical thresholds (see 
text). 



for these kinematic constraints the support is Hmited, as 
Eq.(|5| indicates. A separate consideration, is that inter- 
mediate states carry bare color, so there is no analog of 
unitarity for factorized DVCS. In DIS this distinction is 
irrelevant, but here however, the absence of intermediate 
hadronic states, means that the GPD cannot have the 
proper physical branch cuts. Fig. [T] shows that the gap 
remains even at high . 

The suggestion 5^ that experimental analyses provide 
information only on the kinematical domain on a ridge at 
X — ^ and fixed t and therefore depends on whether 
one can disregard or treat otherwise the unphysical term 
in Eq.(|5| (note that in NLO analyses the domain is 
smeared beyond the ridge [5]). It is this point about the 
sufficiency of the "ridge" that we are examining with care, 
by assuming DRs are satisfied in various model GPDs. 

To illustrate these crucial questions we consider two 
examples of models that should satisfy DRs, namely 
an asymptotic Regge model and a covariant spectator 
model. The Regge pole model contributes to the scat- 
tering amplitude amplitude T(i/, i, Q^) for a single Regge 



FIG. 2: (color online) Threshold induced violations of dis- 
persion relations for both the Regge model (left), and the 
covariant quark diquark model (right), described in the text. 
For the Regge model we show the gap between the calculated 
real parts using the physical threshold at t = O.SGeV^ and 
different values of Q^, and the analytic continuation labeled 
as vth = 0. The same gap is plotted vs. ^ in the lower panel. 
For the quark diquark model calculation the direct (dots), 
the dispersion relation (dashes), and their difference (full) are 
plotted vs. C, for different t values. 

trajectory a(t) in the simple form 

T«(i^, Q^ t) = Pit, g2)(l - e""(*)) (^^) " . (6) 
So the DR should be 

providing that the integral converges. For a low lying 
trajectory or large enough t so that a{t) < 0, this will 
converge. But this relation is exact only for i>th — 0. 
The actual threshold for = is at —t/4M and further 
for non-zero Q^. So the DR is satisfied asymptotically, 
for J/ >> vth- This is illustrated for several cases in 
Fig [2^ where the real and imaginary parts are plotted 
against i' and against (. The directly calculated real 
part and the dispersion relation result for the real part in 
this unsubtracted dispersion relation are quite separated 
for low v or high C. Note that for current typical JLab 
kinematics (Hall B) < 4.5GeV2, \t\ < 2.0GeV2, and 
0.09 < XBj < 0.6, so the non-asymptotic values of 
and t are quite relevant. 

We next consider a quark diquark model with spin- 
less partons (for simplicity, as in Ref. [3 . Because this 
is a covariant model it satisfies the polynomiality con- 
dition thus allowing the GPD to be continued into the 
large ^ or small a;/^ region in which the analyticity re- 
quirements apply [Ij. The subtraction, A(^,t) the dif- 
ference between the evaluation of Eq.([2| and Eq.(|5]) for 
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FIG. 3: Difference between tiie dispersion reiation and direct 
calculation in a scalar quark-diquarfc modei inciuding tlie hard 
jet hadronic mass as from Eqs.( |8|9| . 

the symmetric case is presented in Fig|2]D, which clearly 
displays non negligible ^ and t variations of A(^,i), thus 
demonstrating that A(^,t) cannot be identified with a 
dispersion subtraction constant. In this case, since all 
X, ^ and t dependences are part of the model in a non- 
trivial way, the threshold ^max is necessarily the physical 
one. Given that our subtraction "constant" A{(,t) is ac- 
tually a function of (, due to the threshold dependence, 
we cannot see a direct relation to either the so-called D- 
term or the J=0 fixed pole [8^, although at high \t\ 
and there is a flattening out. We have not addressed 
the nature of the states yet. For GPDs some kind of a 
duality model needs to be introduced that makes the col- 
ored quark and remnant jets look like hadrons (see recent 
study on this subject in addition to providing a pre- 
scription for analytically continuing to the appropriate 
threshold. The prescription we suggest as an alternative 
to analytic continuation aims at reducing the kinemat- 
ical threshold mismatch by replacing the variables used 
in Figjljand Figj2]with variables including a mass, toj, 
for the hard partonic jet. Our prescription is in line with 
Ref. [9] where it was exactly pointed out that kinematical 



threshold mismatches might arise in the collinear fac- 
torization approach if the mass of the hard partonic jet 
is disregarded. Although considering jets with mass is 
not equivalent to hadronization, it might get us closer to 
what a hadronic intermediate state is. Following [21 [TU] 
we replace the hard propagators for the struck quark in 
the hard part of the handbag with a variable jet mass 

1 . , 1 

7 ^ With — 7 . (8) 

The dispersion relation becomes 

r r Hix,(i + '{^)x,t) 

RcH = PV J dXj dm]p{m^j) ^^J^ 

(9) 

where pirri^j) is a jet mass distribution. The results shown 
in Fig|3]demonstrate that the gap obtained as a result of 
having two different thresholds in the massless calcula- 
tion (Figj2| is softened, due to the new set of variables 
that better account for the correct range of integration 
over the partons' virtuality and transverse momentum 
(see also discussion in [3]). 

In conclusion, we have shown the limitations of ap- 
plying DRs to deeply virtual exclusive processes, and 
have given insight into the partonic nature of GPDs by 
examing the role of variables external and internal, re- 
spectively, to the quark loop that appears in the leading 
order factorization formulation. In particular, we show 
that it could lead to misleading results to base global 
parametrizations on DRs as recently done in [T^]. To pin 
down GPDs we advocate comprehensive measurements 
of both the real and imaginary components through var- 
ious asymmetries and cross section components in a wide 
range of all kinematical variables, C, t and Q^. 
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